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TRAJECTORIES WITH THRUST AND DRAG  7\\, 


We abandon now the restriction of zero thrust and zero drag by 
adding air forces and moments and the thrust force. 


Since the air forces depend upon the motion relative to the air 
and since for zero wind velocity the atmosphere rotates with the earth, 
it is advisable to use a coordinate system fixed to the earth. 


In the equations of motion appear mass forces, which are the 
product of the mass and certain accelerations. In the rotating system 
ef the earth there is not only the gravitational acceleration but also 
a centripetal acceleration and a Coriolis acceleration. These additional 
@ccelerations are caused by the so called carrier velocity ( in the case ~ 
at hand the earth is the carrier). They can with vectcr symbolism be 


‘expressed as follows (derivation of the formulae may be found in any 


textbook of analytical dynamics, which uses vectors), 
; Carrier velocity, Uz = [e x -) 
Centripetal acceleration y= (ld x, ) 


Coriolis acceleration a, se bigs (ey xe v ) 


The symbols are defined in Fig. 5, where the orientation of ‘the vectors 
is alse shown. The centri- 
petal acceleration becomes 
zero in the two poles of 
the earth, where the local 
vector and the vector of the 
earth rotation are parallel. 
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Fig. 5. Carrier velocity, v 
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Centripetal Acceleration (4 xv, 
and Coriolis Acceleratior ’ o (Wo ‘ Y) 
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Fig. 6. Forces and Moments on a Ballistic Missile 


=- Drag 


2 
fv 5S Ch (lbs. ) 
ac Sg ORY vs Cy, (lbs. ) 
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Moment about c.g. = @ v* Sifcy CP ibe) 
5) ' 


= Rocket thrust Ci) 

- Weight of missile= m->g (1b. ) 

= Velocity (ft./sec.) 

- Angle of attack (rad. ) 

- Flight Path angle (rad. ) 

= Angle of (Pitch) Attitude = (rad. ) 

= Angle of Control PeMieruton (rad. ) 

= Reference Area for aerodynamic coefficients (ft .<) 
= Hetorenca Length for aerodynamic coefficients ) 
- Mass Density of the air ; (lb « sec*/ft*) 
~ Lever Arm of Control Force Fig ee 3 
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The Coriolis acceleration disappears wherever the missile velocity is 
parallel to the axis of the earth. 


EQUATION OF MOTION FOR PLANE TRAJECTORIES 


We confine our discussions to the most simple case, where the 
motion is contained in a plane. Jt is obvious that this has to be a 
plane which contains the center of the earth. If we want to be exact, 
we have to admit that such a motion can not exist on the rotating 
earth, because the Coriolis force would force the missile out of this 
plane. Let us estimate how great this force is. If one has a velocity 
of 10,000 ft/sec, the Coriolis acceleration can become as high as, 


Cox ¥ = 0.729 » 107% . 10,000 = 0.729 ft/sec* = 0.026 g. 


It is clear that this sideward acceleration cannot be neglected on a 
ballistic missile, where the flight time may be 1,000 sec. and the 
sideward deviation may integrate up to about 20 miles: 


The centrifugal force does not increase with the missile velocity 
and is, therefore, relatively small for a fast missile, The maximum 
centrifugal force occurs at the equator, where one has, 


ee reve 2 g 
Wx(Wxr) = (.729)* ree wig. eh 10 e141 ft/sec* = .0034g. 
This number is an order of magnitude smaller than the Coriolis acceleration. 
We will neglect both accelerations in the following equations. 


The equations of motion are the mathematical expression for Newton's 
law: Mass times acceleration equals the external acting forces. We 
split the forces up into two components, one parallel, the other perpen- 
dicular to the flight path tangent. For the forces parallel to the 
flight path tangent one has 


(37) ™ fi = =—mg ar tha v°S op + Teo I~ aX). 


The term on the left hand side of the equation is easily recognized as 
the acceleration force. On the right hand side comesfirst the component 
of the weight, which tries to decelerate the missile, when the flight 
path angle is positive. Then follows the drag, which also slows the 
missile down and has, therefore, a negative sign. The last term is the 
thrust component in the direction to the flight path. It is assumed that 
the missile is controlled by swiveling the rocket engine or by deflecting 
the direction of the thrust by means of a jet vane. 
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The forces perpendicular to the trajectory are, 


(38) yor. Sn = 
a dt 


aisesa leah + Sc, - Wendi. 


The term on,the left hand side is a centrifugal force caused by the 
curvature, Ct » of the trajectory. On the right. hand side, inside 

the parenthesis, is the term g- saga fl which has the en sign, 
because it tries to make “WME negative. The term v*/(r, +h) is the 
centrifugal force the missile would experience, if it would fly at 
constant altitude around the earth. Fora flat earth, i.e., To =~ 
this term would disappear. If the speed is high enough, this term would 
balance the weight force and one would have orbiting conditions. From 
this condition follows the aren fore for a circular orbit, 


(39) ov Ly g(rot h) Bolo Pe hn)! 


The next term on the right hand side of (38) is the lift, which causes 
a positive 4 /{t and is, therefore, introduced with positive sign. The 
last term is the thrust component, which acts in a direction opposite 
to the lift. 


The balance of moments about the lateral axis through the center of 
gravity of the missile can be written, 


(40) I, Sas ae Stem +t Le dey + Tl,sind 
dt OL 


The term on the left hand side is again the inertia term, the 
product of moments of inertia, Is about the y-axis and the pitch 


acceleration, _& “O. » On the right hand side come first ‘two aerody- 

d : 
namic terms. The first represents the moment one would have, if the flight 
condition would be stationary, that is with zero pitch rate. The second 
term is the aerodynamic pitch damping. Finally there is the moment caused 
by the deflection of the thrust direction, the control moment. The aero- 
dynamic damping term includes the damping caused by the transport of 
propellant masses from the tanks, usually located close to the center of 
gravity, to the outlet of the rocket nozzle. This damping term can be 
written, 


(A)... ye ea ee 
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Where  — is the propellant mass flow per second, e the distance of 


the propellant c.g. from the missile c.g., and s the distance of the 
rocket nozzle outlet from the missile c.g. 


psec Misswir TAT 
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4 Fig. 7. Propellant consumption produces a Pas according 
7 to equation (41). 


Besides the equation of motion one has to consider the control 
equation, (42), which is the equation of motion of the control surface 
or the swiveling rocket engine. This equation may be written, 


a.(7-d, ) + afd-d,) +4, (1-dp) at : 
& an aa 
Ny li-Wp) + Ko(@ Op) Kg lO-G,) Hg (O-G) +h, |@4. ) dd 
lee 


The meaning of this equation has been explained in detail in the lectures 
on missile control. Here it will be sufficient to mention that the left 
hand side is made up of terms which describe the inertia, a5, the damping, 
a1, and the steady state amplification factor, a,, of the control members) 
On the right hand side one finds the different magnitudes, O , 

: i, @ dt, which are fed into the servo, each multiplied by 


"the pokmeupondink gain setting, Ky > K O°? K O° hes’ K,. -» The sub- 


. script p is used to designate the programmed values. 


Equations (37), (38), (40), and (42) have to be integrated in order 
to find the wanted time histories of the variables, 


7G % 4.0%, 


It is obvious that one needs two additional equations, to be able to do 
this. These are the geometrical relation 


® @ = wut 


which can be read from Fig. se and the kinematic relation 


(44) Teale v Sin At 
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REMARKS ON NUMERICAL SOLUTIONS 


The mathematical difficulty of this system of simultaneous 
differential equations is caused by the fact that a number of the 
coefficients of this system are nonlinear experimental functions. 
Such functions are: 


the density of chadyht » & function of altitude, h. 


the drag coefficient, ep, a function of Mach number, angle 
of attack, of , and the Reynolds Number. 


the lift coefficient, cy, a function of angle of attack,  , 
Mach number, aenweet deflection, » if aerodynamic 
control surfaces are used. 


the moment coefficient, c,,, a function of angle of attack, 
o , control deflection, » Mach number, and center 
of gravity position. 


the thrust, a function of altitude. 


If one has to solve the equations of motion for a higher number 
of cases, it is almost impossible to do it without the help of 
electronic computers, Both analogue computer and digital computer 
have the possibility to read arbitrary functions which have been stored 
in them in form of curves or tables or mathematical functions. This is 
a relatively easy procedure, when the functions depend upon one variable 
only. Unfortunately, our arbitrary functions depend upon more than one 
variable. In principle, however, this problem can always be solved by 
curvefitting methods. 


In order to illustrate the conditions a sketch of the curves of the 
drag coefficient are Broueta in 
Fig. 8. The curves are not 

affine so that an assumption 


of the form Ch = F(a) « GiM) 


would not be sufficient. One 


has to use an assumption of» 


the form 
/ Mach No, 
Fig. 8. Drag coefficient, Pe is a (45) or os a | pea) AC) 


function of Mach number ana angle N=/ 
of attack 7 YX, 
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In case of the drag coefficient the series can be discontinued after the 
second or third term. For the moment coefficient, more terms may become 
necessary. 


The various functions needed for the trajectory calculations have 
to be obtained from the Meteorologist, the Aerodynamicist, the Control 
Engineer, and sometimes, when elastic deformations become important, 
even from the Structural Engineer. 


A great number of methods have been developed for performing the 
actual numerical calculations. All methods replace the differential 
equations by difference equation with sufficiently small step intervals. 
Though the use of high speed computers made the time element less 
important than it was, when the computations were performed on desk 
calculators, one is still discussing the question, which is the fastest 
numperical method, The Runge-Kutta method modified by Nystrom and the 
method of central differences are still in use as they were when the 
trajectories for the V-2 were calculated (Ref. 8 and 8a). 
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